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INTRODUCTION 

The dynamic response of large structural systems 
is often analyzed using component mode synthesis (CMS) 
techniques. CMS is widely accepted for predicting 
coupled system response with increased modeling effi- 
ciency and flexibility over conventional methods. CMS 
techniques utilize a reduced set of component modes 
to characterize the overall system behavior. However, 
the inability to adequately model the connections 
between components has limited the application of CMS, 
Connections between structural components, and between 
components and ground are often mechanically complex 
and difficult to accurately model analytically. The 
modeling of these connections can profoundly influence 
predicted system behavior. This is because only the 
connections determine the boundary conditions which are 
imposed upon the system components. Thus, improved 
analytical models for connections are needed to extend 
the applicability of CMS and to improve system dynamic 
predictions. 

Parameter identification ( P I D ) techniques can be 
used to improve predicted response when experimental 
data are available. Modeling accuracy is improved 
with PID by reducing discrepancies between the meas- 
ured characteristics of a physical system with those 
predicted by an analytical model of the system. Many 
techniques are available to carry out this process of 
parameter refinement. Most involve the determination 
of a set of structural parameters which optimally min- 
imize differences between experiment and analytical 
prediction. 

This study explores combining CMS and PID methods 
to improve the analytical modeling of the connections 
in a component mode synthesis model. The approach 
involves modeling components with either finite ele- 
ments or experimental modal data and then joining the 
components with physical connecting elements at their 
interface points. Interface connections in both the 
translational and rotational directions are addressed. 
Once the system model is derived, experimentally meas- 
ured data is used with PID methods to improve the 
characteri zations of the connections between compo- 
nents. Corrections in the connection properties are 
computed in terms of physical parameters. With this 
approach, the physical characteristics of the connec- 
tions can be better understood, in addition to provid- 
ing improved input for the CMS model. 

The identification of connection characteristics 
is simplified by requiring individual components to be 
verified before they are incorporated into the coupled 
system model. This requirement will normally not pre- 
sent any difficulties, since component testing and 
verification has become a regular practice. With this 
requirement, the components are verified before they 
are used in the coupled system model. Any differences 
between the measured and predicted coupled system 


response can be solely attributed to inaccuracies of 
the estimated properties of the connections. Also, 
the quantity of test data that must be obtained from 
the coupled system is greatly reduced. This is par- 
ticularly useful when it is impractical to obtain a 
complete set of vibration test data for a coupled 
structure. Examples, include large space structures, 
spacecraft systems, and turbomachinery. 

Component Coupling Procedure 

Numerous variations of the CMS method are cur- 
rently available for the dynamic analysis of coupled 
structural systems (l to 3). In the classical CMS 
approach, all of the system components are character- 
ized in the modal domain using their respective modal 
parameters (frequencies and mode shapes). Coupling 
between components also is performed in the modal 
domain through use of modal constraints. These con- 
straints are derived from displacement compatibility 
conditions existing at the component interface loca- 
tions. With the classical CMS approach, any compo- 
nents or connections that have been modeled in terms 
of physical coordinates (e.g., finite elements) must 
be transformed into the modal domain before they can 
be included in the coupled system equations of motion. 
The system equations, in terms of modal coordinates, 
are used to compute the system natural frequencies. 

The system mode shapes are computed by transforming 
the mode shapes obtained from the system equations 
back to physical coordinates. 

Recent applications of the CMS method have 
shifted from the classical approach of utilizing only 
modal coordinates. Instead, techniques that use a 
mixture of both modal and physical coordinate systems 
have been implemented (3). There are several reasons 
for the shift to a "mixed" coordinate set. One reason 
is that a combination of component types can be incor- 
porated into the coupled system equations without 
requiring all of the components to be in identical 
coordinate systems. This is particularly useful when 
some of the components have been modeled using F.E. 
methods and other component models have been derived 
from modal test data. In most of the currently used 
CMS methods boundary degrees of freedom of all of the 
components are expressed in terms of physical coordi- 
nates, and the internal degrees of freedom are expres- 
sed in either modal or physical coordinates. The 
inherent efficiency of the component representation is 
retained. With physical boundary coordinates, compo- 
nents can be coupled utilizing classical direct stiff- 
ness assembly techniques as in conventional F.E. 
computer codes. Furthermore, nonlinear connecting 
elements can be used when boundary degrees of freedom 
are in physical coordinates. In the classical CMS 
approach, where modal coordinates are used, it is very 
difficult to incorporate nonlinearities into the cou- 
pled system model because of the difficulties associ- 
ated with defining modal parameters for nonlinear 
elements. 
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This study develops a simplified variation of the 
previously mentioned procedures for CMS. The proce- 
dure is defined to be compatible with PID procedures 
which will be used subsequently for identifying the 
component interface characteristics. The modal compo- 
nents are first converted to "pseudo" finite elements 
to connect modal components to physical finite element 
components. The pseudo elements are then treated in 
the same manner as conventional finite elements, i.e., 
system property matrices are assembled through direct 
stiffness techniques. 

Consider the system shown in Fig. 1. This system 
is comprised of two components which are coupled by a 
physical connecting component. The undamped, free 
vibration equation of motion for the uncoupled system 
is written in terms of physical coordinates as: 
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where [M] and [K] are the component system mass and 
stiffness matrices, and {ul and {u } are the vectors of 
component nodal accelerations and displacements (the 
superscripts refer to the component identification). 
Equation (1) can further be partitioned by separating 
displacements internal to the components from those 
that are at the interfaces between components. When 
this is done Eq. (1) is written as: 
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The coupled system equation is obtained by applying the 
displacement compatibility conditions at the interface 
between the components and the connections. The dis- 
placements of the component and the connection must be 
equal at the interface therefore: 
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r 

i 


u b 

& 


(4) 


From conservation of energy principles and the above 
transformation, the coupled equation of motion is 
found from: 


^coupled ^ ^uncoupled^ 
and 

^coupled = ^ ^uncoupled^ ^ 

Substituting the property matrices from Eq. (2) into 
Eq. (5) the coupled equation of motion is: 
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The procedure outlined above can be used to couple any 
number of physical components into the system equat- 
ions. As mentioned previously, when modal components 
are to be coupled into the system model they are con- 
verted to pseudo physical components and then are 
treated in the same manner as conventional physical 
components. The pseudo physical property matrices are 
obtained from orthogonality relationships between the 
property matrices and the modal parameters. When the 
component mode shapes are normalized so that the modal 
mass matrix equals the identity matrix, the modal and 
physical property matrices are related by: 
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(3) 


U] T [M]M = [i] 


Using Eq. (3), the transformation matrix, [T], which and 

relates the dependent and independent displacement 
sets is: 


[*] T [K][*] = IN& 


(7) 


where [M] and [K] are the component physical mass and 
stiffness matrices, C^w 2 ^] are the component freq- 
uencies, and [<t>] is the matrix of component mode 
shapes. 
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When experimental modal data is used to charac- 
terize the component, the matrix [ <t>] , containing the 
component mode shapes may be rectangular. If "m" mode 
shapes are measured, and the value of the mode shapes 
are recorded at "n" different physical locations on 
the component, then the mode shape matrix will be of 
order n x m. Normally, there will be more measure- 
ment locations available than there will be modes that 
can be measured. To obtain a square modal matrix from 
experimental mode shape data, data at some measurement 
points can be neglected so that the number of points 
is equal to the number of modes. When data at measure- 
ment points is discarded no information is lost as far 
as the overall system response is concerned, so long 
as measurements at the component's interface points 
are retained. Once a square mode shape matrix is 
available, the pseudo physical property matrices are 
related to the modal data by: 

T - 1 - 1 

[M p ] = U T ] M 
and 

[K p ] = [* T ]"Vwl][*] _1 (8) 

where [M p ] and [Kp] are the component pseudo mass and 
stiffness matrices. (The coefficients of the mass and 
stiffness matrices are in terms of physical coordinates 
corresponding to the location and direction where the 
mode shapes are measured). 

The matrices computed in Eq. (8) are designated as 
pseudo matrices because their physical interpretation 
is unlike that of conventional mass and stiffness 
matrices. Because it is impractical to measure all of 
the component modes, the modal data will be incomplete 
(see (4)) and will not contain all the information 
required to produce the actual component mass and 
stiffness matrices. Therefore, although the mass and 
stiffness matrices computed in Eq. (8) are in terms 
of physical rather than modal coordinates, the matri- 
ces will not necessarily represent the actual physical 
mass and stiffness characteristics of the component. 
The mass and stiffness matrices from Eq. (8) will 
reproduce the measured frequencies and mode shapes, 
and will be suitable for representing the component 
in the coupled system model. 

Either the "free" or the "fixed" boundary compo- 
nent modes can be used for the component character- 
ization. The "free" mode shapes are those modes that 
correspond to the component when it is in the uncon- 
strained or free boundary condition. In many situ- 
ations these modes are more conveniently obtained 
than the fixed boundary modes. This is particularly 
true when the modes are measured experimentally, 
because the component itself does not have to be phys- 
ically constrained during the experimental testing. 

In practice, the free boundary condition often is 
approximated by suspending the component from flexible 
cords or by supporting it on soft springs. 

The fixed modes are obtained by simultaneously 
constraining all of the component's boundary degrees 
of freedom while performing the modal testing. Ana- 
lytical! y, the fixed modes are computed as easily 


as the free modes. Experimentally, they are more dif- 
ficult to obtain, because all of the component's 
boundary degrees of freedom must be fully constrained 
during the experiment. To attain this condition 
requires that elaborate fixtures be attached at the 
components boundary locations, and in practice, full 
constraint is never completely achieved. Another dif- 
ficulty of using fixed boundary mode shapes is that an 
additional set of "static" deflection or constraint 
modes must be added to the set of fixed boundary modes. 
These modes are required so that the component will 
have flexibility at its boundary locations where it 
is connected to adjacent components. 

Normally, the values of the experimental mode 
shapes are measured in the translational directions. 

It is not generally practical to measure the values of 
the mode shapes in the rotational directions because 
of limitations in available instrumentation. However, 
it is sometimes desirable to couple rotational degrees 
of freedom between components. If the values of the 
mode shapes are not measured in the rotational direc- 
tions, the pseudo matrices will only have transla- 
tional degrees of freedom and there will not be means 
of coupling the rotational connecting stiffnesses. To 
circumvent this difficulty, the rotational values of 
the mode shapes can be extrapolated from the transla- 
tional values, either by curve fitting through the 
translational degrees of freedom and then computing 
the slope of the curve at the connection location, or 
by us^incj^an approximate F.E. model of the component 

When the rotational values are extrapolated from 
a curve fit any existing rotational inertia effects 
will not be reflected in the values of the rotations. 
Neglecting the actual independent motion of the rota- 
tion implies that there is no rotational inertia and 
that the rotations are dependent on the translations. 
Because of this dependence, the combined translational/ 
rotational mode shapes can not be used directly to 
compute the pseudo matrices without encountering numer- 
ical problems during the matrix inversions in Eq. (8). 

A solution to this difficulty is to initially use 
only the translational mode shapes to compute the 
pseudo matrices. Then, a transformation which is 
based on the dependence between the rotations and 
translations is used to transform the pseudo matrices 
from the translational coordinate system to a com- 
bined translational/rotational system. 

The dependent rotational values of the mode 
shapes can be related to the independent translations 
by: 

«, -i'Ai < 9 > 

J -j 

Where U is the dependent rotation at j, U . are 

0 * A 1 

J 

the translations at the independent measurement points, 
ai are the coefficients relating the independent 
translations to the dependent rotations (determined 
from curve fit, etc.), and n is the number of 
independent measurement points. 

The transformation from the mixed coordinate 
matrices to the entirely translational pseudo property 
matrices is: 
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where [T 1 ] is the transformation matrix derived from 
the relationships in Eq. (9) and {Ua}' is a subset of 
{u*}. For each rotational degree of freedom that is 
added in { u 0 ) * a translational degree of freedom is 
removed from { u A > 1 • The selection of the trans- 
lational degrees of freedom that are removed is arbi- 
trary, and since a translation is removed for each 
rotation that is added, both systems will contain the 
same number of degrees of freedom. 


Using the original translational pseudo property 
matrices from Eq. (8), the transformation in Eq. (10), 
and principles uf conservation of energy, the pseudo 
matrices are derived in the combined translational/ 
rotational coordinate system by: 
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( 11 ) 


Once the component pseudo matrices in Eq. (11) are 
computed, they can be inserted into the system equa- 
tions of motion and coupled to adjacent components 
using the previously discussed procedures. 


where [W] is the weighting matrix and {F} is a vector 
of weighted squared differences. To find the set of 
connection parameters that minimizes the weighted 
squared differences, the derivative of {F} with respect 
to the connection parameters is set to zero. Noting 
that the predicted characteristics (c), are a func- 
tion of the connection parameters (r), the deriva- 
tive of {F} is written as: 

(13) 

Expanding {c} in a Taylor series and truncating higher 
order terms, { c } is approximated as: 

{c} = {c} est + Hrf {Ar} {14) 

Where {Ar} are the differences between the estimated 
and actual values for the connection parameters. 
Substituting Eq. (14) into Eq. (13) and letting 
3 {c}/ 3 {r}=[S] leads to: 

[W] ({c} - { c } e s t - [S]Ur>)[S] = {0} (15) 

From Eq. (15) it is desired to solve for {Ar} so that 
the actual connection parameters can be determined. 
Solving for {Ar} can not be accomplished by simple 
inversions, however, because in general the number of 
measured and predicted characteristics will be greater 
than the number of connection parameters, rendering 
the matrix [S] to be nonsquare. The vector {Ar} can 
be solved for if Eq. (15) is first premultiplied by 
[S]T. When this is done, {Ar} is solved as: 


The final coupled system equations can be used 
to predict the overall system dynamic characteristics. 
The frequencies that are computed from this equation 
will correspond to the overall system resonances. The 
accuracy of the predicted frequencies will be depend- 
ent on the precision with which the connections between 
components have been modeled. It has been assumed 
that the component modal models have been verified 
and are accurate, and also, that the proper component 
modes have been included in the model to adequately 
predict system response (see sample problem one). 

The mode shapes derived from the system equations 
will correspond to the physical degrees of freedom 
included in the system model. When the combined 
translational /rotational model is used some of the 
mode shape values will correspond to translational 
degrees of freedom and some to rotations. The accu- 
racy of the mode shapes, like the frequencies, will 
be dependent on the adequacy of the component modal 
representations and the modeling of the connections. 

Parameter Identification Procedure 

Once the system equations of motion and their 
correspondi ng frequencies and mode shapes are computed, 
and the experimental system modes have been measured, 
PID can be used to find an improved set of connection 
parameters that better predict the measured experimen- 
tal system data. For this study the Weighted Least 
Squares method for parameter estimation is used (6) . 

If {c} and {c} are vectors containing the measured 
and computed system frequencies and mode shapes respec- 
tively, then the weighted squared difference between 
the predicted and measured characteristics is: 

{F} = [W]({C} - <C}) 2 (12) 


{Ar} = ([S] T [W][S]) 1 [S] T [W] ({£} - K } £ST ) (16) 

An updated set of connection parameters is computed by: 

{r} " {r) EST + {Ar} U 7 ) 

or by substituting from Eq. (16): 

{ r} = {r} EST + ([S] T [W][S]) -1 [S] T [W]( { c > - {c} E$T ) 

(18) 

Since {c} is approximated by a truncated series, 
the improved connection parameters will be only an 
approximation to the final parameters. However, the 
final parameters can be obtained by iterating on 
Eq. (18). 

A direct approach for computing the elements of 
the sensitivity matrix [S] is to perturb the analyti- 
cal model with changes in the connection parameters, 
and then compute the resulting changes in the system 
characteristics. The elements are then computed by 
setting Tij equal to the change in the ci charac- 
teristic divided by the change in the rj connection 
parameter. Alternative methods for computing these 
derivatives have been presented (see (_7)) but for 
problems such as the example, with only a small number 
of connection parameters, the above method is adequate. 

The selection of the system characteristics that 
are used in the estimation procedure is determined by 
data acquistion capability. Experimentally, it is 
generally easier to measure frequencies than mode 
shapes, so in many cases it may be practical to include 
more frequencies than modes shapes. Characteristics 
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other than frequencies and mode shapes also can be 
utilized; in (8), it is suggested that kinetic energy 
may be a usefuT characteristic. Once the characteris- 
tics are chosen, the weight that is placed on each 
characteristic must be determined. If one character- 
istic is measured more accurately than another, then 
it can be weighted more heavily. 

When the number of system characteristics is 
large, the size of the weighting and sensitivity mat- 
rices increases, and the matrix in Eq. (18) may become 
ill conditioned for inversion (see (9J). The PID 
procedure only requires a minimum number of system 
characteristics to adequately identify the connection 
parameters since each component has already been veri- 
fied. Therefore, the size of the matrices in Eq. (18) 
will be kept small and inversion problems will be min- 
imized. Another problem may arise when the analytical 
model cannot be exactly made to fit the experimental 
data. When this is the situation the set of connec- 
tion parameters that minimizes the differences, rather 
than eliminates them, must be used. The model may not 
be able to produce the desired measured system charac- 
teristics because of limitations in the component 
modal representation. Also, if the experimentally 
measured modes are not orthogonal, perfect agreement 
can never be achieved because the analytical model can 
only produce orthogonal mode shapes. 

Sample Problem One: Coupled Beams 

The following sample problem is offerred to dem- 
onstrate the component coupling and parameter identi- 
fication procedures. To verify these procedures 
simulated experimental data generated from a F.E. 
model was used. The sample problem (Fig. 2) is com- 
prised of two simply supported beams connected at 
their ends. For simplicity, both beam components 
were made identical. In actual applications the sys- 
tem can be partitioned into any set of components 
that is desired. Each of the components in this 
problem are discretized into seven massless, planar 
beam elements. Concentrated translational masses are 
added between the elements at nodes 2 through 7 and 
10 through 15. The components are connected by a 
rotational spring (K = 10. E5) at nodes eight and nine. 

A connection also is made to ground by a rotational 
spring (K = 10. E5) added to the second component at 
node 16. 

The accuracy of the computed system frequencies 
as a function of the number of modes used for the com- 
ponent representations was evaluated with six, four, 
and two component modes (see Table I). Both the six 
and four component mode representations produced sys- 
tem frequencies that are in good agreement with the 
baseline F.E. solution. Although there are only six 
component modes in the F.E. solution, the six mode 
representation does not produce exact frequencies 
because the F.E. model has more than 6 degrees of 
freedom. The two mode representation allows for the 
first and third modal frequencies to be predicted sat- 
isfactorily but does not provide enough information 
for an accurate prediction of the second and fourth 
frequencies. At least two component modes are required 
so that there will be a rotational degree of freedom 
at each end of the component that is connected to 
ground (only one mode is needed for the other compo- 
nent). In every case the component mode solution 
produced frequencies that are higher than the baseline 
frequencies. This is understandable since the compo- 
nent mode solution uses a truncated set of modes and 
therefore does not include all of the component's 
f lexibi 1 ity. 


For the initial attempt at identifying the con- 
nection properties, only the simulated system frequen- 
cies from the F.E. model (Table II) were used in the 
parameter identification routines. It is preferable 
that the connection properties be identified without 
having to use system mode shapes because the mode 
shapes are considerably more difficult to experimen- 
tally measure than the frequencies. When either six 
or four component mode representations were used two 
possible solutions were found for the K1 and K2 con- 
necting stiffnesses which satisfied the system fre- 
quency constraints (see Table II). The chosen solution 
was dependent on the initial starting estimates for Kl 
and K2. Although neither solution is equal to the 
actual connecting stiffnesses, the first one is rea- 
sonably close considering the limited number of system 
data used and the approximation of the component modal 
representation. When either of the solutions are 
input into the F.E. model they produce system frequen- 
cies that are very close to the exact frequencies. 

The first solution does produce a better set of system 
mode shapes. In an actual application, without more 
than system frequency information, it would be impos- 
sible to determine which of the two solutions is 
closer to the actual values of the connecting stiff- 
nesses. Furthermore, since both the five and two 
system frequency cases produced similar solutions 
there is no advantage to using more than two system 
frequencies. When two component modes are used a 
maximum of four system frequencies are available, 
therefore the five system frequency case cannot be 
analyzed. For the two component modes and two system 
frequency case, the solution failed to converge. 

A subsequent attempt, using a combination of 
both system frequencies and mode shapes was made with 
the expectation that the identification of the connec- 
tion properties would be improved. By adding the 
first mode shape as a constraint, along with the first 
five system frequencies, the second multiple solution 
was eliminated. When only one system frequency and 
one mode shape was used, the problem still converged 
to the first solution regardless of the initial esti- 
mates for the connecting stiffnesses. This combination 
of system data is ideal because, while it eliminates 
the multiple solution, it only requires a minimal 
amount of experimental data. Similar results were 
produced for both the six and four component mode 
representations, while the two mode representation 
continued to present difficulties. 

Sample Problem Two: RSD Rig Verification 

Once the component coupling and parameter identi- 
fication algorithms were evaluated with simulated data 
(Sample Problem One), it was decided to assess the 
procedures using actual experimental data. To accom- 
plish this, the RSD (Rotating Structural Dynamics Rig) 
at NASA Lewis Research Center was selected. The RSD 
rig (Fig. 3) is designed to simulate engine structures 
to study active rotor control and system dynamics 
(component interaction) problems. The rig components, 
although considerably simpler than a real turbine 
engine's, were scaled such that they would simulate an 
actual engine's structural dynamics response 
characteristics. 

The objective of the parameter identification was 
to determine the stiffnesses of the squirrel cage bear- 
ing support that connects each end of the rotor to the 
support frame. To accomplish this, the RSD rig was 
divided in two components; the rotor support frame, 
and the rotor. Each of these components was charac- 
terized verified experimentally, so that accurate 
component representations would be available for the 
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coupled system model. In the system model the support 
frame was represented by an experimentally verified 
F.E. model while the rotor component was represented by 
experimental modal data. Since both components were 
experimentally verified, any differences that appeared 
between the predicted and measured system characteris- 
tics could be attributed to the uncertainties in the 
squirrel cage connections between components. This 
approach considerably simplified the verification task 
by reducing the quantity of modal data required from 
the coupled system. 

The support frame finite element mesh is shown in 
Fig. 4. The frame is mounted on a relatively stiff base 
plate so grid points 35 through 39 are fully constrained. 
Grid points 19 and 20, where the rotor is attached, 
were allowed to freely displace. This free condition 
is representative of the conditions used during the 
modal tests and is also compatible with the require- 
ments for the component coupling procedure. The grid 
points are connected with beam (bending and axial 
deformations) elements except for the diagonal elements 
at grid 35 which are modeled with rod (axial deforma- 
tion only) elements. All of the elements are modeled 
with A36 steel properties. The frame F.E. model was 
analyzed with NASTRAN, to compute the component fre- 
quencies and mode shapes (Fig. 5). The frequencies 
were experimentally verified by using vibration data 
obtained from an HP 5423 Dynamic Analyzer. The rotor 
modal representation was obtained by measuring the 
rotor mode shapes in the free boundary condition. This 
condition was approximated by hanging the rotor from 
bungy cords. The component modal characteristics were 
generated from transfer function data obtained from the 
dynamic analyzer and impact testing. A total of six 
rotor modes were measured (see Fig. 6} including two 
rigid body and four elastic modes. 

The support frame and rotor were coupled by com- 
bining the physical F.E. model of the frame with the 
modal representation of the rotor. For simplicity the 
coupled system model was constrained to motion only in 
the vertical plane. This restriction allowed for a 
reduction in the required number of degrees of freedom 
in the system model and allowed for all of the system 
testing to be performed in one plane. The coupled sys- 
tem frequencies for the six mode rotor representation 
are plotted along with the measured frequencies in 
Fig. 7. The predicted frequencies were computed for 
different values of squirrel cage stiffness to deter- 
mine the effect that the cages have on the system fre- 
quencies. To generate these results it was assumed that 
both squirrel cages had identical stiffnesses. This was 
a rational assumption, since both cages are built to the 
same specifications. (Subsequent to this analysis the 
cage stiffness was measured as 5050 lb/in. using a 
static loading test. ) 

Only the first three computed system frequencies 
are shown because only three frequencies were measured. 
When all three frequencies are used the cage stiffness 
is identified as 5750 Ib/in. This value is in good 
agreement with the measured stiffness (5050), consider- 
ing that only three system frequencies were used for 
the parameter identif ication. In Fig. 7 it is shown 
that this amount of difference in cage stiffness does 
not have a significant effect on the system frequencies. 

In addition to the six mode rotor representation, 
a four and two mode representation were used to deter- 
mine the effect that the number of component modes has 
on the stiffness identification. The four mode repre- 
sentation identified the same cage stiffness as the six 


mode representation. The two mode representation iden- 
tified the cage stiffness as about 2300 lb/in. or only 
46 percent of the measured stiffness. It was expected 
that the two mode representation would be insufficient 
for identifying the cage stiffness because this repre- 
sentation is inadequate for accurately predicting the 
system modes. It is obvious that the two mode repre- 
sentation cannot produce very good results because 
only rigid body modes are included in the representa- 
tion, and the system modes involve elastic bending in 
the rotor. Although rules of thumb are available for 
determining the required number of modes, additional 
work is required in this area. 

CONCLUSION 

From the two sample problems analyzed in this study 
it was determined that the stiffness characteristics of 
component connections can be identified using component 
mode synthesis and parameter identification procedures. 
Furthermore, the characteristics can be identified using 
experimentally obtained component modal representation 
and a minimal quantity of measured system modal data. 

In the first sample problem it was found that multiple 
solutions are possible, but that they can be avoided 
when system mode shapes are included in the identifi- 
cation procedure. In the second problem it was found 
that the rotor for a rotor/support frame coupled sys- 
tem could be adequately represented by experimentally 
obtained modal data. It was also found that only three 
system frequencies had to be measured for the connec- 
tion characteristics between the frame and rotor to be 
identified. From the results obtained thus far, it is 
determined that the quantity of data required for the 
component representations and for the connection char- 
acteristic identification is problem dependent. There- 
fore, each application must be treated on an individual 
basis. 
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TABLE I. - COUPLED SYSTEM FREQUENCIES (SAMPLE PROBLEM ONE) 


(Connection stiffness k-| = k2 * 1 Oxl 0 5 ] 
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finite- 

element 

solution, 

Hz 



Number of component modes 
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Component mode synthesis solution 



Eigenvalue 

Frequency, 

Eigenvalue 

F requency , 

Eigenvalue 

Frequency, 


+ percent 

Hz 

+ percent 

Hz 

4- percent 

Hz 

4.94xl0 5 ' (112) 

5.11x10 s 

*2 

114 

5.1 2x1 0 4 
*2 

114 

5.25x10 s 

4-3 
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1 . 07x10 s (165) 

1.23xl0 6 

♦7 

177 

1 .28xl0 6 
4-9 

180 

1 .57x10 s 
4-21 

197 

6.99xl0 6 (421) 

7 . 1 3x1 0 fe 
+1 

425 

7 . 1 7x1 0 6 
4-1 

426 

7.54x10 s 

4-4 

431 

9 .72x1 0 s (496) 

1 .08xl0 7 
+1 

523 

1 . 09x1 0 7 
4-1 

525 

1 .72xl0 7 

561 

3.39xl0 7 (927) 

3.43xl0 7 

932 

3.44X10 7 

933 
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TABLE II. - COMPUTED CONNECTION STIFFNESS 


Number of 
system 
frequencies 

Number of component modes 

6 

4 

2 

Connection stiffness 

k l 

k 2 

k l 

k 2 k 1 k 2 

5 

6.5x10 s 

7.4xlC 5 

6.3x10 s 

7.3x10 s 

(a) 

3.7x10 s 

13.1x10 s 

3.6x10 s 

12.5x10 s 

2 

6.7x10 s 

7.4x10 s 

4.1x10 s 

10.8x10 s 

(b) 

3.7xT 0 s 

13.4x10 s 

5.4x10 s 

8.1x10 s 


a 0nly four system frequencies available. 
^Solution does not converge. 
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FIGURE 2.- COUPLED SYSTEM (SAMPLE PROBLEM ONE). 
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MODE 1: 0 Hz 



MODE 2: 0 Hz 



MODE 3: 141 Hz 



MODE 4: 304 Hz 



MODE 5: 609 Hz 
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SQUIRREL CAGE STIFFNESS, lb/in. 

FIGURE 7.- COUPLED FRAME/ROTOR ANALYSIS (6 
COMPONENT MODES). 
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